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A MATHEMATICAL MODEL OF THE PHENOMENON OF RADIOACTIVE FALLOUT 


R. R. Rapp 


To describe the phenomenon of fallout we must, of course, begin with 


the nuclear detonation. The main distinguishing feature of a detonation . | 


of this type is a tremendous release of energy almost instantaneously and 


within @ very small space. The actual amount of energy in calories is 


1 


related to kilotons of yield in the ratio ‘0! to 1 according to the 


Atomict Energy Commission book Effects of Atomic Weapons. Or, in other 


words, one kiloton is equivalent to one trillion calories of energy. ’ 


To get a meteorological comparison for. this amount of energy, we 


might compute the energy released by the condensation of water. Suppose 


we consider a depth of rainfall of one centimeter over one square kilo- 


meter, or, 10+° gms of water. Allowing 5hO calories per gram of water 


12 calories, or approxi- 


condensed, the energy released vould be 5.4 x 10 
mately 5 kilotons. The big difference, of course, is the fact that the 
reinfall releases its energy throughout a large volume of atmosphere and 
over a period of time soacired in minutes or hours, vhereas the bomb | 
releases its energy within a very limited space and in a matter of seconds. 

This release of energy creates a hot bubble of air which rises because ere 
of its buoyancy. The cheracteristic mushroom-ehaped cloud associated with 
nuclear explosions is formed by water vapor condensed as a esniie of 
vertical lifting. The temperature gradient between the hot bubble and the | 
surrounding atmosphere causes & strong toroidal circulation. 


Another unique feature of a nuclear explosion is the formation of 


fission fragments. These are the end products of the splitting of the 


753 992 


‘Uranium atoms and are almost all radioactive. Fission fragments, together 


vaith any other debris mixed into the fireball at an early stage, are carried 
aloft, and aa the cloud rises, some of the material may be left behind 
‘forming @ wake or stem. The cloud will stop rising when the energy 
lavailable to lift it has dissipated. At this time, which 1s referred to 

as the time of stabilization, the spatial distribution of the radioactivity 
provides one set of the initial conditions required in the fallout model. . 
| Within this stabilized cloud are vitually all the fission products , 
‘formed by the nuclear detonation. Radioactive particles may be formed by 
‘condensation of vaporized material or by the coalescence of fission products 
with solid or liquid earth particles. The size of the particles is deter— 
‘mined largely by the initial conditions of the Surat 3} if little solid | 
material is taken into the cloud, the particles will in all likelihood be 
emall; if, on the other hand, the device ia detonated on the surface, 

many tons of earth will be swept into the cloud and a sizeable fraction 

of the fission products will be contained in large particles. The size 

of the particles containing radioactive elements has an important bearing 
won the fallout preblem. The radioactivity located in or on large particles 

will br brought rapidly to the ground, that lodged in small particles will 

remain suspended in the atmosphere for long periods of time. Another of 

the initial conditions required in the fallout model is, therefore, the 

distribution of activity with particle size in the stabilized cloud. 

One of the features of the fission products which mist be considered | 

is the rate of decay. If tha: eats of decay is such that little radio~ 
activity is left when the particles reach the ground, the effect on life 


will be considerably lessened. The elements decay according to an expon~- . 


ential law, each with its own characteristic half life. When the decay. 
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rates of all the different elements are added together, the decay of the 


mixture is given approximately by a power law. 


“1.2 | 
A(t)= A(1)t_ | a 


/ 
Because of this decay it is necessary to carefully consider the time when $1 
‘particles reach the ground in computing the pattern of dose rate: 
Figure 1 shovs the distribution of activity with particle size as 
reconstructed from two fallout patterns guesived after surface bursts . 
“Glididw ie Pancave? “uneaveatnty GP aacadinqaent, “ule Alatesbution-Gppears 
‘to be valid for both small devices in Nevada and large devices in the 
Pacific. The limited amount of data available did not provide the means 
foe estimating the effect of soil type on particle size distribution. 
Given the stabilized cloud and the associated distribution of activity 
“with particle size and in space, it is necessary to determine the position 
of the activity after it falls to the ground. The rate of fall is of 
primary importance here. The first analytical method which came to mind 
for describing rate of fall was Stoke's Lav, but it was soon found by 
examining fallout samples that some large particles wae in the range of 
one thousand microns diameter, and Stoke's Law does not apply to such 
large particles. ‘Therefore, an aerodynamic law of fall was used which 
takes into account, by means of drag coefficients, inertial as well as 
viscous forces. Figure 2 shows the comparison of different fall rates | 
based on experimental data. The solid curves are based on Stoke's Law 
for 10, 50-and 500 micron radius particles. Note that Stoke's Law is 


not sensitive to the density to the atmosphere; thus the variation of 
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FIG.I— RADIOACTIVITY AS A FUNCTION OF PARTICLE SIZE 
FOR THE STEM AND MUSHROOM OF AND ATOMIC CLOUD 
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HE ATMOSPHERE 


AS A FUNCTION OF HEIGHT FOR VARIOUS PARTICLE SIZES 
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| fall velocity with height is due ‘solely to changes in viscosity. The ' 
dashed curves » representing the serodynamic law for perfectly smooth | 
spheres, show that the fall rate is sensitive to the density of the atmos— 
phere and further that various particle sizes are affected quite differently. 
Nov, the actual particles involved are by no means smooth spheres; thus, it 
is not quite proper, although we have done this in our calculations, to 
use the smooth—sphere law. The dotted curves show the fall velocities based 
on drag coefficients measured and reported in the Handbuch der Seperinen 
Physic for what the experimenters term "riindlicher Korper." 

Of eourse, as the particles start to fall through the atmcephere, they 
will be carried by the wind. There are a grest many wind effects which have 
to be taken into aceount. First, the whole cloud will be moved bodily by 
the mean wind at the cloud level. At the same time, there will he shearing 
effects due to the variation of the wind with altitude and turtmlent effects, 
small scale eddies, which will tend to spread and tear the cloud. Figure 3 
shows the wind effect. This picture, prepared by Dr. Jamea Edinger of UCLA, 
shows cubic clouds in order to demonstrate best how the different sizes of . 
eddies affect the clouds. Sketch (a) shows the simple translation due to 
the mean wind. Sketch (b) shows what would happen if there were eddies of | 
the same diziension as the aloud, preserving it fairly well, but turning 
and twisting it. Bketeh; (c) ‘shows that if there were eddies of an order of 
magnitude exiller than me ‘cloud, little wisps of debris would be pulled 
out and wisps. of clear aig. avates sais The effect of. shear is show in 
sketch (a) - ae cloud will be arevn ato @ leng ribbon. If all these 


motions of the atmosphere are put togetber, sketch (@) shows what the aquare 
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(A) Gross movement— atmospheric motions of a (B) Horizontal deformation— atmospheric motions 
scale large compared with cloud dimensions of a scale comparable to cloud dimension 


OY 


Height 


45,000 ft 


to 35:000Ft aoe) 


(D) Vertical deformation —shear of 
horizontal wind in the vertical 


(C) Eddy diffusion— atmospheric motions of a 
scale small compared with cloud dimensions 


Horizontal! section 


Vertical section 


(E) Total effect of gross movement, 
deformation,and eddy diffusion 


Fig. 3 


‘cloud would look Like after a period of time under the combined influence 


a all these motions — ~ not exactly 4 pleasant thing to try to deal with 


‘mathematically. 


; 


It must be obvious without going into the detaile that there are many 


‘facets of this problem which are not well understood. The distribution of 


‘activity with size is not well known. The spatial distribution of the 
‘activity has not been messured, and the small—scale motions of the atmos— 
phere present a statistical fluctuation which can not be predicted. ioueven: 
in order to make an approximation, & simplified picture has been constructed 
‘to try to explain fallout. Figure 4 show the cloud in model form. It 
rises to height H,; H, is @ fimetion of the amount of energy released and 
the ambient atmospheric conditions at the time. It 1s assumed to be a 


‘eylinder (the mushroom) of @ certain diameter and thickness set on a stem | 
of about 1/5 this diameter. 

As far as the activity distribution within this cloud is concerned, it 
tas been assumed that most of the activity, because of the great stability 
of the smoke-ring circulation, will rise into the upper part, with a small 
amount trailing down into the stem. It is assumed that there is a constant 
mixing ratio through the mushroom. This means that 4 count of particles 
per unit volume will show an exponential decrease through the mushroom. It 
has been assumed also that the concentration is constant throughotit’ the 
stem. ' 

As for the distribution of activity with different particle sizes Ae 
recourse to some classified results has ylelded a rough idea of the activity 
on particles of various sizes. This is shown in Fig. 1. For fall velocities 


@ choice of laws exists. The roundish body results reported in the literature 
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FIG.4A—ASSUMED DISTRIBUTION OF 
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_geem most consistent with the physical picture, but to simplify the 
| , : ; 
‘calculation, the smooth sphere results of Pig. 2 were used. 


. With these simplifications it is possible to make some arithmatic 
( 


, Computations and produce & calculated fallout sahara which will serve 


ae ee 


: for & comparison with observed results. In The Effects of Atomic Weapons): 
: Appendix FP, the equations for the MeBecvorres of particles are shorn. 


' These have been modified slightly and put into the following form: 


Qs) \ ; 


- Where x is the distance out along some coordinate, r is the radius of 


particle, H is the height at which it starts, Wis the fall velocity and 


~ 


V 1a the wind velocity. Note that the turbulent diffusion has been 
‘ignored; it 16 assumed that the effect is small compared with the effects 
‘ef shear and translation. We hope that diffusion will, in part, counter-— 
balance errors in our assumption of @ cylindrical cloud. 


The cloud wes divided into increments of height, and each slice vas 


divided into increments of particle size. The trajectories for each of 
these groups of particles wes then computed for the center of the cloud. 
This ebiadline has been programmed for the IBM-70l1, which provided for a 
‘ hundred different intervals of height and a hundred different intervals of 
_ particle size. The position formulae can be applied to get the xy position 
of the center of any given group of particles of size r from height H; 


from the activity distribution, the radioactive Pission product associated 


with these particles is known. The groups of particles are agcumed to 
’ maintain a circular pattern. When all the positions on the ground are 


. known, the contributions from all of the different (rH) groups to the 
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| radioactivity at a series of points can be summed. The pattern of radio~ ' 
t 
{ 


: 


“active fallout can be drawn from the data at these points.‘ 

In order to show the tremendous effect of the wind on the patterns of 
! fallout, Figs. 5 and 6 are presented. These two patterns were caloulated ; 
| using the model just deseribed. The parameters which are functions of .: 
| bomb yield have been held constant, and the winds have been changed. 


iy 
Figure 5 shows the patterns at 1/2, 1, 2, and 6 hours for a moderately 


, strong wind with little shear. Figure 6 shows the patterns for the same 

‘ timeg but under conditions of light and variable winds. ‘These patterns - 

, ave probably extreme cases, but they were both calculated from actual 

wind soundings and serve to demonstrate the critical nature of the 

_ meteorological factors of the fallout problem. | 

| The model is crude, and it would not be reasonable te expect that by 

| assuming certain facts about the size, shape, and distribution in the 

| cloud an accurate prediction could be made for an actual erent: However, | 

7 patterns so constructed are similar to those that are found after test shots. 

‘This leads us to believe that we are on the right track physically and that 
we will be able to improve our model and knowledge of fallout by following 

' along this track and trying to get the necessary data to replace our 


assumptions with, observed facts. 
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